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Abstract: Infinite order differential equations have come to play an increasingly signif- 
icant role in theoretical physics. Field theories with infinitely many derivatives are ubiq- 
uitous in string field theory and have attracted interest recently also from cosmologists. 
Crucial to any application is a firm understanding of the mathematical structure of infinite 
order partial differential equations. In our previous work we developed a formalism to 
study the initial value problem for linear infinite order equations with constant coefficients. 
Our approach relied on the use of a contour integral representation for the functions under 
consideration. In many applications, including the study of cosmological perturbations in 
nonlocal inflation, one must solve linearized partial differential equations about some time- 
dependent background. This typically leads to variable coefficient equations, in which case 
the contour integral methods employed previously become inappropriate. In this paper we 
develop the theory of a particular class of linear infinite order partial differential equations 
with variable coefficients. Our formalism is particularly well suited to the types of equa- 
tions that arise in nonlocal cosmological perturbation theory. As an example to illustrate 
our formalism we compute the leading corrections to the scalar field perturbations in p-adic 
inflation and show explicitly that these are small on large scales. 
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1. Introduction 

Applications of infinite order differential equations to theoretical physics have attracted 
considerable interest recently. Such equations are ubiquitous in string field theory (l| (see 
Q for a review) and also arise in a number of toy models of string theory such as the 
p-adic string H [|| and discrete world-sheet models ||. Moreover, such equations have 
recently attracted interest from cosmologists due to a wide array of novel cosmological 



behaviours ||-[22]. Of particular interest are nonlocal inflationary models [23]-[^7]], such 
as p-adic inflation, which can provide a unique playground for studying string cosmology in 
an ultra-violet complete setting. Nonlocal inflation has the remarkable property that slow 
roll inflation can proceed even when the potential is naively too steep and may therefore 
offer a novel way to circumvent the difficulty of finding flat scalar field potentials in string 
theory. This remarkable behaviour was first observed in [23] and subsequently generalized 



to a broader class of models in |26| (see also The perturbative analysis of 23]-[26| 

was verified using fully nonlinear simulations in p7| . 

Perhaps the most interesting feature of nonlocal inflation is the possibility of generating 
an observably large nongaussian contribution to the temperature anisotropies in the cosmic 



microwave background. In [24] an estimate was provided for the nongaussianity in a wide 



class of nonlocal models and this estimate was verified by a more quantitative calculation 
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in [pql . It is worth noting that the prediction of [24] is consistent with a subsequent 
claimed detection of nongaussianity in |28|| . The quantitative calculation of [25] was made 
possible by a special gauge choice and by working to only to zeroth order in the e slow 
roll parameter, despite keeping up to first order in the r\ parameter. This is a consistent 
approach for the models under consideration due to the hierarchy \rj\ ^> e, which is typical 
in hill-top inflation models [29|. However, more generally one would like to be able to 



study cosmological perturbations in nonlocal inflation without making such assumptions. 
Progress in this direction is stymied by the fact that the nonlocal cosmological perturbation 



theory equations are extremely difficult to solve. In [3C] the theory of linear infinite order 
equations with constant coefficients was developed. 1 However, this theory relied on integral 
transform methods which fail to be useful for the variable coefficient equations that arise 
in cosmological perturbation theory. 

To improve on the calculation of p5| , then, clearly we require more sophisticated 
analytical tools and a better understanding of the formal aspects of infinite order equations. 
In this paper we make progress in this direction by developing the theory of a broad class 
of infinite order variable coefficient equations. Our formalism is particularly well suited to 
studying the kinds of equations that describe nonlocal cosmological perturbations. This 
work is the first in a series of papers that aim to develop a general approach to nonlocal 
cosmological perturbation theory. As an illustration of our method, we will compute the 
leading corrections to the inflaton perturbations in p-adic inflation, showing explicitly that 
these are small and do not lead to super-horizon evolution. 

It is appropriate, before moving on, to mention the difficulties and complications that 
arise when one wishes to interpret higher derivative theories as fundamental. Such theories 
are often fraught with classical instabilities known as Ostrogradski instabilities |37]] (see 
|38| ] for a modern review), as emphasized in |39[j . In this paper our primary interest is in de- 
veloping the tools necessary to solve nonlocal equations of motion in a very general context 
and we will not address the important question of when such theories can be phenomeno- 



logically viable. See [22, 40 1 for examples of stable, interacting nonlocal theories. See [41] 
for a discussion of the nonlinear stability of theories with infinitely many derivatives. 

The organization of this paper is as follows. In section || we introduce the class of 
equations under consideration and review previous results. In section || we consider the 
simplest possible extension of the analysis of [3C]. In section || we give a formal solution 



of the class of equations under consideration which relies on knowledge of the inverse of 
the nonlocal operator under consideration. In section || we give two different methods for 
computing this inverse operator. In section [6| we apply our methods to p-adic inflation. Fi- 
nally, in section |7| we conclude by discussing further possible applications of our formalism. 
In appendices A and B we review the solutions for the inflationary background and scalar 
field perturbations in p-adic inflation. In appendix C we discuss some technical manipula- 
tions involving Bessel functions. In appendix D we discuss the non-self-adjointness of the 
d'Alembertian in de Sitter space. 



1 See also different approaches to the initial value problem, [|35| for mathematical analysis of 

p-adic and string field equations and |3g] for more details on solving nonlocal equations. 
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2. Infinite Order Equations with Variable Coefficients 



In this paper we develop the theory of a particular class of linear infinite order partial 
differential equations having variable coefficients. The class of equation under consideration 
can be cast in the form 

F(D)0(t,x) = m 2 (^(i,x) (2.1) 

where □ = g^V^Vj, is the covariant d'Alembertian. 2 Throughout we will assume that the 
kinetic function F(z) (also called the generatrix in the mathematics literature) is entire. 
With this assumption F{z) can be represented by a power series centered at z = and 
having infinite radius of convergence. Thus, we define the action of the pseudo-differential 
operator F(p) on some smooth function <j> by the series expansion 

oo 

F(D)</>=^a n D> (2.2) 

?1=0 

where the powers are, of course, understood as composition of the operator □ with 
itself n times and the coefficients in the expansion are 

a n = : — (2.3) 

n! 

Our restriction to analytic F(z) is merely for simplicity. We expect that it should be 
straightforward to generalize our results to F{z) having isolated poles (such as the zeta 



strings model [|42j1 ) or branch cuts (see [30] for mathematical analysis and |43] for cosmo- 
logical applications). 

The motivation to consider ( |2.lD comes from studying a general class of infinite order 
equations 

F(n)<f> = V , [<f>] (2.4) 

which are typical in string field theory and also in cosmological models. Suppose one wishes 
to solve ( |2.4p for small inhomogeneities 8<j)(t, x) about some known homogeneous solution 
<po(t). Plugging the ansatz 

0(t,x) = 6,(t)+£0(t s x) (2.5) 



into (^J) and linearizing in 6(f) we obtain 

F(p)6<i>(t, x) = V" IMt)} 5<Kt, x) (2.6) 



which is precisely of the form Q2,I| ) with m 2 (t) = V" [(f)o(t)]. 

Equations of the form ( |2.I| ) with m 2 = const and g^ v = if*" were considered in our 
previous work pJ. In that paper solutions were derived using the formal operator calculus 



and we were able to exhaustively count the number of initial data necessary to uniquely 



specify a solution. The formalism developed in [3C] relied on the fact that in flat space 



with m = const (|2.]|) is an infinite order equation with constant coefficients. However, in 



2 Our choice of metric signature is such that g^ v = rf v = diag( — 1, +1, +1, +1) in flat space. 
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many applications both g^ v and m 2 will depend nontrivially on space-time coordinates. In 
this case the contour integral approach adopted in (3(1] is no longer applicable. 



A physical example where an equation of the form (2.1) arises is the computation of 
cosmological perturbations in nonlocal inflationary models. In the case of p-adic inflation 
the dynamical equation for the scalar field is 



P 



-□/(2m2) 



(2.7) 



For cosmological applications one should solve yQ. 



ds 2 = g^dx^dx" = -dt 2 



in an FRW geometry 
- a 2 (t)dxidx l 



Slowly rolling background solutions (f>o(t) have been constructed which source a quasi-de 
Sitter expansion H = a/a = const in |23|]. Neglecting induced inhomogeneities of the 
metric, the cosmological perturbations satisfy the equation 



-0/(2mJ) 



P 



<ty(t,x) = P 4>l~\t)-l <ty(*,x) 



which is precisely of the form (2.1) with 

F(U) = 
m 2 (t) = 



D/(2m 2 B ) . 
(*) 



(2.9) 



(2.10) 
(2.11) 



3. The Case m 2 (t) = 

Let us first consider the simplest possible extension of the results of |3(| by assuming that 
m 2 (t) = but allowing for g^ u to be nontrivial. 3 Let us assume that the function F(z) 
has iV zeroes, all of which are order unity (this is to be expected for physically interesting 
operators). 

To solve (^^) we first provisionally assume that <fi is a formal eigenfunction of □, 
meaning that it is a non-trivial solution of the equation 

U(f)=-L0 2 (j) (3.1) 

Acting on (j) with the full operator F(n) we have 

oo 

H u )<f> = 

n=0 

oo 

= J> n (-„ 2 )> 

n=0 

= F(-lo 2 )(/) 

3 This analysis applies equally well to m 2 (t) — const since any additive constant could be absorbed into 
the definition of F(D). 
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(The re-summation on the third line is justified for any value of uj because F{z) is entire.) 
Hence the solutions of (|3.1| ) will give rise to solutions of ( |2,1[ ) provided the eigenvalue co 2 
is chosen to satisfy the transcendental equation 

F(-uj 2 ) = (3.2) 

By assumption this equation has N distinct roots u> 2 (i = 1, • • • ,N) and hence there are 
iV distinct solutions of (2.1). We assume that to each of these roots, there corresponds 



a formal eigenfunction fa of the d'Alembertian. 4 The most general solution of Q2.1| ) is 
obtained by superposing these eigenfunctions as 



N 



</>(t,x) =Y J Mt,*) (3-3) 
i=i 

Each fa is a solution of the second order equation Ufa = —ujjfa and hence contains two 
degrees of freedom. 5 The full solution ( |3.3|) then admits 2N initial data, two for each zero 



of F{z). This is precisely the same result that was obtained for flat space in [30 where 
a physical interpretation was provided in terms of the poles of the propagator. It should 
be remarked that we are not making any completeness assumption on the set of formal 
eigenfunctions. In fact, the construction being local, we need not even assume at this stage 
that the wave operator □ is essentially self-adjoint. 

Let us briefly discuss the application of this approach to p-adic inflation. Consider 
equation (|2.9|). During inflation fa(t) = 1 so that m 2 (t) <C 1 and ( EM ) can be approximated 
by 

6fat,x)=i0 (3.4) 



p-o/«_p 



which is of the form under consideration. This equation is solved, as above, by taking 
<5</>(i,x) to be an eigenfunction of □. This is precisely the approach that was adopted in 
|23f| to study cosmological perturbations in p-adic inflation (and also in [p4| , in a more 
general context). In this paper we will to go beyond this approximation and show how to 
systematically include the effect of having fa(t) different from unity. 



4. The Method of Successive Substitution 

Let us now develop the theory of equation ( |2.1[ ) in the case where the time dependence 
of m 2 (t) cannot be neglected. 6 Let us suppose first that we are able to solve the simpler 
equation 

F(D) x = (4.1) 

4 There is of course a scaling ambiguity in the choice of 4>i- In the context of the cosmological applications 
considered in this paper, the ambiguity will be resolved by considering solutions with prescribed asymptotics 
at infinity and at the horizon (that is simply the usual Bunch-Davies procedure, see Appendix B). 

5 By "degrees of freedom" here we refer to the freedom to specify two independent functions of the spatial 
variables x, to specify the Cauchy data corresponding to a solution. This is equivalent to saying that, upon 
quantization, each 0, will have its own set of annihilation/creation operators. 

6 In fact, the approach of this section may fail in the case of constant m 2 (t). 
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This falls into the class of equations considered in section || and the solutions are eigen- 
functions of □. Denoting the N solutions of (4.1) by Xi we have 



□Xi = -"iXi (4-2) 
F(-cof) = (4.3) 

for i = l,--- ,iV. 

Next, consider the inhomogeneous equation 

F(p) x = J (4.4) 

Let us suppose that we are able to solve this equation to obtain the particular solution 
Xpar- 7 Writing the particular solution as 

Xpar = GJ (4.5) 

defines the operator G, which is the inverse of the kinetic function F(0). In the mathe- 
matical literature G is often referred to as the resolvent generatrix. 

We now demonstrate that knowledge of the solutions of ( |4.1[ ) and of the inverse operator 
G is sufficient to construct the full solutions of (|2.1|). Consider the infinite series 



= Xi 

+ G [m 2 Xt ] + G [m 2 G [m 2 X i\] + G [m 2 G [m 2 G [m 2 X i\]} 

+ ■■■ (4.6) 



Plugging (|4.6D into fl2.1| ) and making use of the fact that F(0) Xi = and that G is the 
inverse of F one can easily verify that 4>i affords a formal solution of Q2.1|). The most 



general solution of (2.1) is obtained by superposing the modes 4>i as 

N 

^(t,x)=J^^(t J x) (4.7) 

i=l 

Hence we conclude that <f>(t,x) contains the same number of degrees of freedom as x(t, x )- 
That is, (f>(t,x) admits two initial data for each zero of F(z). 

The utility of the solution ( |4.6[ ) lies in the observation that when m 2 (t) is a small 
perturbation the series of successive iterations on the second line takes the form of an 
expansion in powers of a small parameter. 8 In this case the series (|4.6| ) can be truncated 
at some finite order to obtain an approximate expression for 0j(i,x). 

Notice that the definition of the quantities F(z) and m 2 (t) in ( |2.1| ) is ambiguous. The 
form of this equation and also the expression for the solution (4.6) is unchanged under the 



substitution F(z) — > F(z) + A, m 2 (t) — * m 2 (t) + A for any constant A. In applications, 
we can take advantage of this ambiguity to define m 2 (t) in such a way that it is small and 
the series ( |4.6|) can safely be truncated at low order. Indeed, we have already implicitly 



7 We will discuss the solution of this equation in more detail in the next section. 
8 We will quantify what is meant by "small" on a case-by-case basis. 
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done this in writing d2,9|) : by subtracting p8(j) from both sides of the equation we have 
defined F(n) : m 2 (t) in such a way that m 2 (t) <C 1 during inflation. With this definition 
the application of (|4.6[ ) to equation ( |2.9[ ) has a nice physical interpretation. The solution 
Xi represent the pure de Sitter space modes constructed in [23] while the subsequent terms 



in the series (2.9) represent the leading order corrections that arise due to the slow motion 



of (j>o(t) away from the unstable maximum <f>= 1. 



5. The Inverse Operator 



In the last section we showed that ( |4.6| ) affords a formal solution of (2.1). This solution 
is given in terms of the solutions of the simpler equation ( f4.1|) and the inverse operator, 
G. In section || we showed that the construction of the solutions of is straightforward 
once the eigenfunctions of □ are known. Hence, the success of our method hinges on our 
ability to construct the inverse operator G associated with F(p). In |3(| we were able to 
obtain exact expressions of the solutions of the inhomogeneous equation for the special 
case g^ v = r)^ v by using contour integral methods. However, generalizing these results to 
curved backgrounds is a nontrivial task. Below we discuss two approaches to computing 
the resolvent which will prove to be useful in applications. It is worth emphasizing that 
the methods developed here are applicable to a broad class of inhomogeneous nonlocal 
equations. 

5.1 Expansion of the Source in Eigenfunctions 

Let us first consider the case where G acts on an eigenfunction of the d'Alembertian. Hence, 
we wish to solve (0) when J satisfies 

Uj = m 2 jJ (5.1) 



If we assume that F(m,j) ^ then (fOj) is trivially solved by 

Xpar = / 2 \ J ( 5 ' 2 ) 

r (rrij) 

It follows that, when G acts on an eigenfunction of □, we have the simple expression 

G(a) = — * (5.3) 
F(m z j) 

Suppose, now, that J is not an eigenfunction but it can be expanded into a sum of 
eigenfunctions as 

J = ^c n J n (5.4) 

n 

where the c n are constant and □</„ = m^J n . Note that we are not assuming that the J n 
form a complete set of eigenfunctions of □. All we require at this stage is that J should 
be expressible as a series in a set of formal or true eigenfunctions of □, which should be 



-7- 



differentiable term by term infinitely many times. Assuming that F(m^) ^ for all n we 
have 

owing to the linearity of the resolvent G. We will see that the simple expression ( |5.5| ) will 
allow us to compute the corrections to the inflaton perturbations in p-adic inflation. 

Before moving on let us comment on the generality of this approach. The form (|5.4| ) 
is quite special and we are not, in general, guaranteed that an arbitrary source J(t, x) is 
expansible in eigenfunctions of □. This is so because for nontrivial geometries the operator 
□ may fail to be self-adjoint when one imposes physically interesting boundary conditions 
on the eigenfunctions. For example, in de Sitter space the d'Alembertian is not self-adjoint 
when acting on a function space that includes Bunch-Davies normalized mode functions; 
see appendix D. Of course, even when □ is not self-adjoint it may still happen that some 
particular physically interesting source term is expansible in eigenfunctions, which is the 
case for p-adic inflation. 

5.2 The Method of Infinite "Differentiation" 

As we have discussed, the case considered in the last subsection is rather special and an 
expansion of J into eigenfunctions of □ may not always be possible. Thus, it may be 
interesting to develop another approach to computing the resolvent which applies more 
generally. Here we develop an alternative formalism which is based on a modification of 
the "method of infinite differentiation" employed by Davis [ft4f| . The method developed in 
this subsection will not actually be required for the explicit example we consider below and 
is included for completeness. 

The method of infinite "differentiation" will furnish us with an expression for £?(□) as 
an infinite series expansion in powers of □ . The price for the generality of our construction 
will be that the series cannot in general be summed to obtain a closed form expression for 
G. For reasons that will become clear shortly we assume that F(0) ^ (so that an 7^ 0). 9 

Now, let us proceed with the construction of the resolvent. Our aim is to solve the 
inhomogeneous equation fl4.4p . To this end, consider the infinite series of equations formed 
by operating on ( |4.4| ) with successively higher powers of the operator □. That is, consider 
the series of equations 



F(n) x = J 
nF(p) x = uj 
a 2 F(a) x = a 2 j 



9 This does not entail any loss of generality since if F(0) = then we can redefine F(z) — » A + F(z) and 
2 — > A + m 2 (where A is any constant) without changing the form of (2.]). 
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Writing F{z) in terms of the expansion ( j2.2| ) we then have the system of equations 

a x + aid* + a 2 a 2 X + a 3 n 3 x + • • • = J 
+ a a x + aiD 2 x + a 2 a 3 x +•••=□ J 
+0 + a D 2 x + ai Q3 X + ••• = ° 2 J 



One can now consider (5.6) as an infinite system of equations in infinitely many un- 
knowns X) a X> d 2 X> " ' ' This system can be solved algebraically for the unknown x- The 
equations ( |5.6| ) can be considered as a matrix equation 



AV = S 



(5.7) 



where the matrix A is 



A 



while V is the vector of unknowns 



a 


a i 


fl'2 


«3 





a 


ax 


a 2 








Of) 


ax 











«o 



(5. 



r 



a 2 x 
a 3 x 



(5.9) 



and the vector S is constructed from the source 



S 



J 

□ J 
□ 2 J 
□ 3 J 



(5.10) 



Notice that the matrix A is triangular (this is so because the coefficients a n defined by 
Q2.3| ) are constants). This allows for a straightforward iterative inversion of (|5.6|). Let us 
now illustrate how this works. As a first step we use the first line of ( |5.6| ) to solve for x m 
terms of J and □ n x with n > 1: 

X= — [J-a 1 Dx-a 2 n 2 X + ---] 
a 



Next we use the second line of (|5.6|) to eliminate in favour of □ J and n n x with n > 2, 
giving 



X 



1 

« 



J 



o,i / a 

— aj + 



Of) 



1 
a 



a 2 n z x + 
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Continuing in this matter we obtain an expression for <j) solely in terms of □ n J: 



X 



1 

a 



J - — □ J + -o (a? - a 2 a )n 2 J + 



Writing \ = GJ we obtain a formal expression for the resolvent (?(□) as a power series 
expansion 

oo 

G(n)=J2Kn n (5.11) 



n=0 



where the explicit expressions for the first few coefficients are 



-4(a? ~ a a 2 ) 



1 



— (2a a 1 a 2 -a(- a a 3 ) 



(5.12) 
(5.13) 

(5.14) 

(5.15) 



At the formal level, this procedure can be continued to arbitrarily high order in The 
validity of the procedure will of course depend on the choice of source term J. This is in 
principle a potential limitation of the method. In fact, there are instances in which the 
expansion of x i n powers of □ J fails to converge pointwise. 



6. Application to p-adic Inflation 



6.1 Setting Up the Calculation 

We now apply the formalism of sections |] and || to study the dynamics of scalar field 
perturbations during p-adic inflation, equations (|2~9- 2.11 ). The homogeneous background 
solutions (f>o(t), H(t) were constructed in [23] and these result are reviewed in appendix 



A. The solutions are written as a series expansion in powers of the small parameter u(t) 
defined by 

u{t) = a(t)M = e^ Hot < 1 (6.1) 



where 



2m| 
3Hl 



(6.2) 



is a slow roll parameter satisfying |^| <C l 10 and Hq is the background Hubble scale (written 
explicitly in terms of model parameters in appendix A) . To leading order the quantity m 2 (t) 
(defined explicitly in (2.11)) is given by 



m 2 (t) = -p(p-l)u(t) 



(6.3) 



10 In p-adic inflation the spectral index is given by n s — 1 = 2r/. For the WMAP5-preferred value n s = 0.96 
m we have 77 = —0.02. We will use this value in our examples below. 
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Because m 2 oc u it follows that the correction terms on the second line of fl4.6|) are controlled 
by the smallness of u. If we work exclusively to leading order in u we need only consider 
the first term on the second line of (|4.6| ), 

In what follows it will be simplest to work in terms of conformal time r defined by 



adr = dt in terms of cosmic time t. To leading order we have a = e 



H t 



-1/(H t) so 



that u(t) = 

(-Hot) - '"!. To the same accuracy the d'Alembertian takes the form 



□ 



-T 2 d 2 T + 2rd T + r 2 V 2 



(6.4) 



To proceed with the construction of the solution (|4.6[) we need two ingredients: the 



functions Xi defined by ( [4.l|) and an expression for the inverse operator G, defined by ([4.4) . 



Let us first discuss the solutions of QO|). The relevant functions Xi were derived in |23| 
and these results are reviewed in appendix B. We consider only the tachyonic excitation 
with effective mass — Uq = —2m 2 . and thus drop the subscript i on the solution. Expanding 
X in terms of annihilation/creation operators ak, a k and c-number valued modes functions 
Xk we have 

d 3 k 



X(*,x) 



akXfc(^)e ik ' x + h.c. 



(6.5) 



(2tt) 3 / 2 

where h.c. denotes the Hermitian conjugate of the preceding term. To leading order in u 
and Tj the mode functions can be written as 



Xk(r) 



(-krf'H^ (-kr) 



(6.6) 



up to an irrelevant constant phase. For ease of notation we introduce the Fourier space 
d'Alembertian operator defined by 

□ fe = -H 2 [r 2 d 2 - 2rd T + (kr) 2 ] (6.7) 

By construction the function fl6.6| ) obeys the eigenvalue equation 

O k Xk(T) = -2m 2 x k (r) (6.8) 

Our goal is now to compute the corrections coming from the second line of ([4,6|), To 
leading order in u equation ( [4.6[) takes the form 

Mr)-Xk(r) = G[J k (r)} (6.9) 

where we have introduced the notation 

Jfc(r) = m 2 (T)xk(T) 

= -*-i)l^(f)''(- fcT » 3/2+ "<-,(- fcT ) < 6 - 10 ' 

Obviously to proceed further we will need an expression for the inverse operator G. We 
could proceed by expanding G as a power series in □, following ( 5.11| ), and then simply 
computing □?Jj.(t) to high order in n. Unfortunately, this will turn out to give rather 
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poor convergence properties for the resulting series (this observation does not imply that 



G [Jfc(r)] is divergent, only that the approach of subsection ^2 is not ideal for this particular 



problem). The failure of the formalism from subsection 5.2 in this physically interesting 
case is disappointing, however, the simple trick discussed in subsection ^ will be quite 
sufficient for our purposes. 

6.2 The Large Scale Limit 

As a warm-up exercise we first consider the limit — kr — > which is relevant for cosmological 
observations. In this limit the "source" term J^(t) given by eqn. ( |6.10| ) obeys the eigenvalue 
equation 



□ fc J fc (r) -Ami ( 1 + ) Jk{r) 



(6.11) 



To prove this notice that in the limit — kr — > the function x becomes (see appendix B 
for more details) 

H 



Xk(r) 



so that, from ( 6.10 ), we have 

Jk(r) 



-kr) 



p(p-l)-^=(kH G T{-Tfv 



(6.12) 



(6.13) 



Now, notice that the operator — > —Hq(t 2 8 2 — 2rd T ) in the same limit. Acting with 
this operator on fl6,13|) we have U k J k {r) = —&\t]\Hq{\ + 2 1 77 1 /3 ) Jfc(r). Using the definition 
of the slow roll parameter, r\ = —2m 2 / '(JSHq), we trivially recover the result ( |6.11 ). 

Given the result ( |6.1lD and the method of subsection 5A it is easy to compute the 
correction term on the right-hand-side of (|6.9|): 



G[J k (r)) 



1 



F 



-4m 2 . ( 1 + 



2|*7l 



Mr) 



1 



2(l+2|r,|/3) 



-Jk(r) 



P 



(6.14) 



(See equation 



Thus, we conclude that the solutions of (2.9) take the form 



Mr) = Xk(r) 



p l+4|r,|/3 _ 1 



u{t)+0(u 2 {t)) 



(6.15) 



in the limit — kr — > 0. We see that on large scales and to leading order in u the slow motion 
of (f)Q yields a tiny almost-constant correction to the solutions obtained in |2^| . Significantly, 
we see that slow motion of the background does not lead to any super-horizon evolution 
for the inflaton perturbations (at least at the linearized level and to first order in u), as 
one would expect in the absence of entropy perturbations [fT^ |. 

6.3 Corrections to the p-adic Mode Functions 



Now we compute the right-hand-side of ( |6.9D without taking the large-scale limit. To 
do so, we first decompose the "source" term Jfe(r) - eqn. ( |6.10| ) - into a discrete sum of 
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eigenfunctions of O k . n The technical details are discussed in appendix C. The result is 
(see equation ( p-8| )) 

Jk(r) 



-/>(/>-!)( -j- 



E^ 1) f^(-^) 3/2 ^/ 2+2 .(-^) 



n=0 



^^(-kr)^J_ 3/2+2v+2n (-kr) 



(6.16) 

n=0 

where the coefficients a$ depend only on r\ and are explicitly defined by ( C-9| , |C-10| ). Each 
member of the sum ( |6.16| ) is an eigenfunction of as can be seen explicitly from equation 
( C-3| ). To compute the right-hand-side of ( |6.9| ) we use formula (5J3). The result is again a 



sum of the form ( 5.16 ), only the coefficients have changed: 



G[J k {r)} = -p{p-l){-^ 



,n=0 



+ E^ 2) f^(-^) 3/2 ^3/2 +2 , +2 .(-^) 



n=0 



(6.17) 



(See appendix C for more details, 
defined explicitly in (|C-lgjC-19|) . 



The coefficients (3n \ @n in the expansion ( |6.17 ) are 



We now plug the result ( 5.17 ) into (|6.9| ) and cast the leading order mode functions in 
the form 

Mr) = Xk(r) [l + u(r)A( 



-kr)} 



(6.18) 



In ( 6.18 ) we have defined the function A(— kr) = G [Jk{ T )} / i u ( T )Xk( T )] which represents 
the coefficient of the 0{u) contribution to fractional difference between (/>fc(r) and Xfc( r )- I n 
figure [I] we plot the modulus | A(x)| on large scales x < 1. For illustration we take rj = —0.02 
and p = 5. In this figure we have retained up to n = 25 in the summation ( |6.17 ), however, 
keeping only up to n = 1 the curve would have been almost indistinguishable. We have 
verified both analytically and numerically that 



lim A(x) 



P 



1 



p l+4|»,|/3 _ l 



(6.19) 



in agreement with ( |6.15 ). 

In figure [2] we plot |A(x)| for a larger range of x in order to show the large scale 
behavior of the correction. Again we take rj = —0.02 and p = 5 for illustration. Here we 
see that at large x the function A(x) undergoes oscillations of amplitude p. Notice that 
even though u <C 1 it may still happen that up = 0(1) before the end of inflation, provided 
p is sufficiently large. In that case our result ( |6. 18| ) for </>fc(r) cases to be reliable on very 
small scales because higher order terms (such as the second and third terms on the second 
line of (|4.6|)) become important. This should not be taken as evidence that 4>k( T ) differs 



11 As mentioned previously and discussed in appendix D, such an expansion is not, in general, possible. 
Our success here relies on the particular form of the source term defined by (|6.1C|) . 
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significantly from Xk( T ) on small scales when p > 1. It is easy to see that this is not so. 
For time intervals At <C Hq 1 one can treat u(t) as a constant. Denoting this constant 
value by u c we see from ( |2.9| ) that <p obeys the equation 

[F{a k )+p(p-l) Uc ] 4> k (r)=0 

The solutions 4>k( T ) are given by eigenfunctions of where the eigenvalues are given 
by the zeroes of F(—lo 2 ) + pip — l)u c . Recall that Xfc( T ) is an eigenfunction of with 
eigenvalue given by the zeroes of F{—uj 2 ). We conclude that on small scales the solutions 
0fc(r) of (|2.9D differs from Xfc( T ) only by a tiny correction to the effective mass. 



1.25 



1.20 



1.15 



1.10 



1.05 



1.00 



0.95 




Figure 1: |A(ir)| versus x on small scales 
x = — kr < 1. Parameters were chosen as 
rj = —0.02, p = 5 for illustration. 



Figure 2: |A(x)| versus x, showing large 
scale behaviour. Parameters were chosen as 
r] = —0.02, p = 5 for illustration. 



7. Conclusions and Future Directions 



In this paper we have developed the theory of a class of variable coefficient equations of 
infinite order which arise when studying perturbations about a time-dependent background. 
We have shown that solutions can be efficiently computed in terms of the eigenfunctions 
of the d'Alembertian and the inverse operator associated with the kinetic function F(0). 
We have illustrated our method for the case of p-adic inflation by computing the leading 



order corrections to the inflaton perturbations derived in [23]. We have shown explicitly 
that these are small and do not lead to any super-horizon evolution. 

It is worth emphasizing that many of the techniques which we have developed here 
- in particular the method of successive substitution in section |] and the construction of 
the nonlocal Green function in section || - are quite general. We believe that even ignoring 
the main themes of this paper such techniques constitute a valuable contribution to the 
understanding of infinite order differential equations. 
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We conclude by discussing some further possible applications of our formalism. One 
motivation for this work was to develop the tools necessary for a fully rigorous and system- 
atic approach to nonlocal cosmological perturbation theory. In the most general scenario 
the solution of the perturbed scalar field equation will proceed very much analogously to the 
calculation in section || only the de Sitter space d'Alembertian is replaced by a perturbed 
d'Alembertian which includes also metric inhomogeneities. This replacement significantly 
complicates the construction of the eigenfunctions and also of the resolvent generatrix. 
However, no new conceptual obstacle is involved since our formalism does not make any 
specific assumptions about the detailed form of the metric. The correlation functions of 
the curvature perturbation may then be computed using the Seery et al. formalism for 
working directly with the field equations E^j. 

Of course, the applications of our formalism are not limited to nonlocal inflation. 
This approach could also be used to study perturbations about rolling tachyon solutions 
in string field theory, such as the Hellerman and Schnabl |4j| solution. We expect also 
that our approach could be straightforwardly generalized to study nonlocal dynamics in 
theories where F{z) has poles or branch cuts. Yet another potential application of our 
method is to the construction of more general nonlocal inflationary background solutions, 
for example large field inflation models which might give rise to observable gravitational 
wave signatures. We intend to return to these applications in future work. 
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APPENDIX A: Review of p-adic Inflation Background Solutions 



In this appendix we review the slowly rolling solutions obtained in |23| for p-adic inflation. 
Consider p-adic string theory coupled to Einstein-Hilbert gravity: 



where 



In 



S 



9p 



I 



Ml 



-□/(2m2), 



+ 



1 



P+l 



m s 



(o/) 1 / 2 is the string mass and we have defined 

2P- 1 



9 2 P = //: 



p- 



(A-l) 



(A-2) 



(A-3) 



with g s the open string coupling constant. The Lagrangian ( |A-2|) is derived for p a prime 
number, however, the theory can analytically continued to any integer value. 
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In [23] inflationary solutions (fto(t), a(t) of the theory ( A-l ) we constructed by employ- 
ing an expansion in powers of u = e A< 



4>o(t) = \-^<j) r e 

r=l 
oo 

H(t) = H - H r e r 



(A-4) 
(A-5) 



r=l 



where, of course, H = a/a = dta/a. We have chosen the parametrization of the solutions 
such that at t — > — oo the field (ft starts from the unstable maximum of its potential, (ft = 1, 
and the universe undergoes a de Sitter expansion with Hubble constant Hq. As t increases, 
the corrections terms e rXt become more important and the field rolls towards the true 
vacuum (ft = 0. The origin of time is chosen arbitrarily which, of course, has no impact on 
any physical observable. 

The solutions obtained in 123], up to order e 2A ', are 



1 



H = Ho 



e \v\H t + 



1 



p — 1 



2 p l+2|r?|/3 _ l 
i _ plnp p-hl e2MHot 



2 p- 1 

In writing ( [A-6 - A-7 ) we have defined the background Hubble scale 

m\ p — 1 



(A-6) 
(A-7) 



m 



6M2 9 2 (p + 1) 



and the slow roll parameter 



2mj 



(A-8) 



(A-9) 



which satisfies |?7| ^ 1. Notice that the small parameter u(t) can be written in terms of 
the scale factor a(t) as 



u(t) 



MHot 



a(t) 



(A-10) 



To leading order in u we can consistently treat the background expansion as pure 
de Sitter space while still working to nontrivial order in the r\ slow roll parameter. The 
consistency of this approach is due to the fact that there is a large hierarchy between the 
77 slow roll parameter (eqn. |A-9| ) which controls the size of the inflaton (effective) mass 
and the e parameter (e = —H/H 2 ) which controls the departures of the geometry from de 
Sitter. In the same approximation we can write the quantity m 2 (t) ( |2.11| ) for the p-adic 
perturbation equation (|2.9|) as 

m 2 (t) -p{p - l)u(i) = -p(p - l)e mot (A-ll) 
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APPENDIX B: Construction of x f° r p-adic Inflation 



In this appendix we consider the solutions of equation ( |4.1| ) for the case of p-adic infla- 
tion. In this case the kinetic function is given by ( 2.10| ) and hence the equation under 
consideration is 



-□/(2m2) 



x = o 



fB-1) 



p 1 v sJ — p 

This equation belongs to the class considered in section ||. The solutions Xn of (|B-1| ) are 
eigenfunctions of □ 

□ Xn = -UnXn 

The eigenvalues solve the transcendental equation F(— ui\ 



(B-2) 

and are given explicitly by 



2 ± 



4irin 
lnp 



(B-3) 



-2m2. 



with n = 0, 1, • • • The n = mode is the usual tachyon with effective mass — oj 2 
The infinity of n > states have complex mass-squared and are ghost-like (contributing 



negative kinetic energy to the Hamiltonian) ; see |27[],|ll[] and |50|]. These states are, pre- 
sumably, artifacts that would not be present in the full string theory. Here we simply omit 
the ghost modes and focus our attention on the tachyon. Since we only consider n = we 
will drop the subscript n on the eigenfunctions. 



The n = solution of (B-2) in a de Sitter geometry is well known 

d 3 k 



(2vr) 3 /2 



«ke ikx Xk(i) + h.c. 



Xk(t) = y 



7T 



a 3 H 



A: 



aH 



(B-4) 
(B-5) 



where aj~, at are annihilation/creation operators, Xk are c-number valued mode functions, 
h.c. denotes the Hermitian conjugate of the preceding term and a = e Hot is the scale factor. 
In ( |B-5| ) the real-valued constant 5 is an irrelevant phase the order of the Hankel functions 
is 




(B-6) 



to leading order in the rj slow roll parameter. 

It is convenient to introduce conformal time r related to cosmic time t as adr = dt 
In terms of conformal time the scale factor is 

1 



a(r) 



H t 



so that r runs from — oo to as a goes from to +oo. Notice that k/(aHo) 
small parameter u in ( |A-1C| ) can be written as 



u{t) = (-H r) 



(B-7) 
hr. The 

(B-8) 
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Small (sub-horizon) scales corresponds to — kr S> 1. By construction, the solutions 
( |B-5| ) have the following small scale asymptotics 

Xk — — 7= for — kr — > oo (B-9) 



a V2A; 

This normalization corresponds to the usual Bunch-Davies vacuum choice. Large (super- 
horizon) scales corresponds to — kr <C 1. The solutions (|B-5| ) have the following large scale 
asymptotics 

H 



Xk 



(—kr) 11 for — /cr — > 



(B-10) 



The factor A; 3//2 corresponds to the exactly scale invariant part while the factor (—kr) 11 
give the slight departure from scale invariance with spectral index n s = 1 + 2rj. 



APPENDIX C: Eigenfunction Decomposition of Jfe(r) 
In this appendix we would like to show that the quantity Jkij)-, defined by 

J k (x) = A k (-krf/ 2 ^H^ 2 _ v (-kr) 

where 



(C-1) 



(C-2) 



(see eqn. 6.10) can be expanded into a discrete sum of eigenfunctions of the operator 



in de Sitter space, defined by (6.7). 

Before proceeding, let us briefly review some facts about the eigenfunctions of It 
is straightforward to verify the identity 



-kTf' 2 C a {-kT) 



I -a 2 ) II 



o 



{-kTfl 2 C a {-ki 



(C-3) 



which is valid for C a any solution of Bessel's equation with order a. In particular, equation 
Q is valid for C a = H { ^,H^,J a , Y a (or any linear combination thereof). 

As a first step to decompose Jfe(r) into a sum of eigenfunctions of then, we should 
write x^ilg^^^) as sum of terms of the form C a (x) (here we have denoted x = —kr). 
To this end we write the Hankel function in terms of Bessel functions of the first kind as 

x^H^/ 2 _ v (x) = x* [J 3/2 - v (x)+iY 3/2 ^(x)] 

= (l + iC0t[(3/2-77)7T]) x"J3 /2 „„(x) 



3/2+jj 



(*) 



sin [(3/2 - 77)71-] 

(this follows from the definitions of Hv and Y v ). Next, we employ the identity 

— v 00 
(l) M X ) = ^c{n\ (l,v)J,j, + 2n{x) 

c{n;n,v) 



71=0 



H + 2n T(n + n)T(u - \i + 1) 



n\ r(y - (i + l- n)Y(v + l+n) 



(C-4) 

(C-5) 
(C-6) 
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which is derived on page 139 of |4q| . Using ( |C-5[) we can we-write the quantities x v J 3 / 2 _ ri (x) 
and x 11 J-3/2+t){ x ) appearing on the second and third lines of ( |C-4j ) in terms of Bessel 
functions of the first kind. The result is: 



3/2+V H 3/2- v ( x ) = 2" [1 + i cot [(3/2 - r,) tt]] £ c(n; 3/2, 3/2 - rj) x 3 / 2 J 3/2+2n ( 3 



3/2-tj 
2»» 



n=0 



" sin [(3/2 -qH ? C(n; " 3/2 + 2?? ' " 3/2 + ^ X3/2 J ~3/ W2n(*) (C-7) 

Using the identity (C-7) we can finally decompose ( |C-1[ ) into a sum of eigenfunctions of 
as follows: 



Hr^ r > 



(1) 



.n=0 



2kV 2 



(-krf 2 J 3/2+2n {-kr) 



n=0 



where the coefficients depend only on n: 



«l 2) ^(-^) 3 / 2 J_ 3/2W2 n(-^) 



(C-8) 



Q W = 2" (1 + < cot [(3/2 - rj) tt]) c(n; 3/2, 3/2 - rj) 
2 V i 

c(n; -3/2 + 2r), -3/2 + rj) 



a 



(2) 



(C-9) 
(C-10) 



sin [(3/2 - rj) tt] 

Equation QC-S| ) is the main result of this appendix. The terms in the summation ( 6.1G| ) 
should be compared to the de Sitter modes (|6.6|). This shows that each term in the 
expansion ( |C-8| ) behaves very much like the mode function for a massive field in de Sitter 
space. 

If we introduce the notation 



(i) 



-Hq-v/tt , 
2k 3 / 2 



(-krf 2 J 3/2+2n (-kr) 



(C-ll) 
(C-12) 



for the eigenfunctions appearing in the summation ( |C-8| ) then, from ( |C-3| ) , we have the 
identities 



V4 2) = ^(-krf 2 J^ /2+2v+2n (-kr) 



□ fc ^) = [-4n 2 -6n]H 2 4n } 
□ fc ^ 2) = [-4(n + r?) 2 + 6(n + ??)] tf 2 



Using equation ( |5.3[ ) we can easily compute G ip^ (r) 



G 
G 



1 



(2) 



F [(-in 2 - 6n) ^] 
1 



. We find 



F[(-4(n + r ? ) 2 + 6(n + 7 ? ))i/ 2 



(2) 



21 



(C-13) 
(C-14) 



(C-15) 
(C-16) 
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where F(z) is defined by ( |2.10| ), Since the resolvent G is linear we have 



G[J k (r)] = -p(p-l) 



V k 



n=0 



(-fer) 3 / 2 J3/ 2+2n (-A:r) 



+XX 2) |gf(-^) 3/2 ^3/ 2Wn (-^) 

n=0 



where the coefficients are 



(i) 



j,3|tj| l j 



a 



(2) 



p5 j^[2(n+u) a -3(rH^)] _ p 



In deriving (|C-18| , |C-19D we have used the fact that r\ = — 2m 2 /(3ifo). 



(C-17) 



(C-18) 
(C-19) 



APPENDIX D: On the Non-Self- Adjointness of the d'Alembertian 

In this appendix we review some basic facts about Sturm-Liouville theory to demonstrate 
that the operator □ fails to be self-adjoint in a de Sitter geometry when the function space 
contains the standard Bunch-Davies mode functions. We are interested in the eigenvalue 
equation 

□v TO (r,x) = m 2 f m (r,x) (D-l) 



In a de Sitter geometry where □ is given by ( |6.4| ). The equation for the Fourier modes can 
be cast in the standard Sturm-Liouville form as 







L [v k>m (r)) + m 2 a 4 (r) v k>m (r) 
where the Sturm-Liouville operator is 

L = d T [a 2 (T)8 T ] +k 2 a 2 {r) 
the eigenvalue is A = m? and the weight is a 4 (r) (with a(r) = — 1/(H t)). 



(D-2) 



(D-3) 



Although ( D-3 ) is in Sturm-Liouville form, it may fail to be self-adjoint if we impose 
physically interesting boundary conditions on the eigenfunctions. To see explicitly how this 
failure occurs, let us suppose that our function space contains the functions f k (r), g k (r) 
defined by: 



where 



fk(r) 
9k{r) 

v f = i 



krf 2 H^{-kr) 
{-krf 2 H^(-kr) 



4"W 



'9 mj 



(D-4) 
(D-5) 

(D-6) 



- 20 - 



The functions /fc(r), g k (r) correspond to de Sitter space mode functions normalized ac- 
cording to the Bunch-Davies prescription. For simplicity we assume that < m/,m s < H, 
however, our analysis does not rely on this assumption in any crucial way. 

Working on the interval r G (ti,T2) the standard theory gives the identity 

T2 drt k {T)L [g k (r)} - H drg k (r)L [f k (rf] 



[a\r) (f k (rrd T g k (r) - g k (r)d T ^(r))]^ (D-7) 



Let us now show explicitly that term on the second line of ( p-7| ) does not vanish when 
one takes the physically sensible endpoints t\ = — oo, T2 = 0. Using the known large- 
and small-scale asymptotics of the de Sitter space mode functions (see appendix B) it is 
straightforward to show that 

[« 2 (r) (r k (r)d T g k (r) - g k (r)d T tf(r))] - -i (D-8) 

in the limit r — > — 00. On the other hand, we have 

[a 2 (r) (A.(r)^ r5fc (r) - 5fc (r)c\/ fc »)] - (-fcr)-"/-- (D-9) 

in the limit — r — ► 0. Clearly ( p-8|) does not equal ( p-9[) and this completes our proof that 
L is not self-adjoint. 

It is a simple exercise to see that, because 

° drr k (r)L[g k (T)}^ f drg k (r)L [f* k (r)} 
-00 J — 00 

the usual theorems about orthogonality and completeness of the eigenfunctions Vk jTn (r) 
on the interval — 00 < r < do not hold. Hence, there is no reason to expect that a 
general source J k (r) in the function space will be expansible in a series of eigenfunctions. 
Of course, this analysis does not rule out the possibility that some physically interesting 
source terms, such as ( |6.10| ), may happen to be so expansible (see appendix C). However, 
this is a special property of the particular source defined by ( |6.10D , rather than a generic 
property of the eigenfunctions of □ in de Sitter space. 

Note that physically the failure of L to be self-adjoint is easy to understand. It occurs 
because the mode functions have very different behaviour on large scale as compared to the 
small scale asymptotics. Hence we expect similar behaviour in general FRW space-times. 

It is interesting - if tangential to our main line of inquiry - to consider more carefully 
the reason that □ is not self-adjoint in de Sitter space. The choice of Bunch-Davies mode 
functions is closely tied to the choice of conformal coordinates defined by 

ds 2 = a 2 (r) [-dr 2 + dx l dx l ] (D-10) 

with — 00 < r < and a(r) = —1/(Hqt). But it is well-known that these coordinates do 



not cover the whole de Sitter space-time [51]. Let us reconsider this exercise using global 



coordinates [51] defined by 

ds 2 = -dP + 4cosh 2 (t)dx^Xi (D-ll) 
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Taking — oo < t < +00 these coordinates cover the full de Sitter space-time. One could 
reconsider the solutions of the eigenvalue equation □% = m 2 v m in this coordinate system, 
however, this effort is not necessary. In the asymptotic regions t — > ±00 the metric ( p-11 ) 
corresponds to a FRW universe with scale factor a(t) ~ e ±Hot . Hence, we can normalize 
the eigenfunctions such that 

- C xl^H^ ( 4r) ( D -i 2 ) 



2 V a 3 H u \aH , 

asymptotically as t — > ±00 (see eqn. B-5 ). Due to the symmetry under t — > —i we have 
same kind of behaviour for the modes in each asymptotic region (that is freeze-out). Now 
there is no inconsistency with imposing periodic boundary conditions on the eigenfunctions 
at t\ = —T and ti = +T. These boundary conditions render □ self-adjoint and one can 
send T —* 00 at the end of the calculation. 



However, the possibility of using (D-ll) (or, indeed, any global covering) to render 
□ self-adjoint is quite irrelevant for our purposes. Certainly the use of global coordinates 
is inappropriate for the physical problem at hand. First off, the covering (D-ll) contains 
an unphysical contracting region for t < 0. Moreover, we are not actually interested in 
real de Sitter space-time but rather inflationary spaces that merely mimic de Sitter for 
some number of e-foldings. Hence, for the physically reasonable choice of coordinates and 
vacuum the non-self-adjointness of □ is simply something we need to live with. 

The reader who still finds the non-self-adjointness of □ disturbing may wish to consider 
an analogy. Using the standard Minkowski coordinates ds 2 = —dt 2 + dxidx 1 one has, in 
Fourier space, oscillatory eigenfunctions e~ rujt . We can ensure that □ is self-adjoint by 
imposing periodic boundary conditions at ±T and sending T — > 00 at the end of the 
calculation. But now consider the same exercise using Milne coordinates. The Milne 
universe looks like an expanding FRW space-time with scale factor that grows linearly in 
cosmic time, however, it is quite equivalent to ordinary flat Minkowski space-time after a 
coordinate transformation. In conformal coordinates the Milne metric takes the form 

ds 2 = a 2 (r) [-dr 2 + dr 2 + sinh 2 r (d8 2 + sin 2 6dc/) 2 )] (D-13) 

with — 00 < r < 00 and a(r) = e T . These coordinates do not cover the full space-time 
but only the patch t 2 — Xi% % > 0. The analogue of the Bunch-Davies modes ( p-4 , ET5| ) are 
solutions of the form 

/ A (r) = Jl—±-— J_ lX ( m/e -) (D-14) 

V * y smh(7rA) 

9x(r) = a/| , . =. = J-iX (m g e T ) (D-15) 

V * y sinh(7rA) 

corresponding to the conformal vacuum |]52|]. (Here the real number A is the eigenvalue 
of the spatial Laplacian and plays the role of k in the de Sitter modes ( P-4 , D^5| ).) The 
modes ( D-14 , p^l5 ) behave very differently near the endpoints r = ±00 and we encounter 
the same difficulty that we had in the de Sitter case. The boundary conditions implied by 
the conformal vacuum are inconsistent with the self-adjointness of the d'Alembertian in 
Minkowski space. 
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